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It  has  recently  been  established  that  the  local  spin  density  approximation 
(LSDA)  within  the  density  functional  theory  can  provide  an  accurate  theoretical 
estimate  of  static  structural  properties,  phonon  spectra,  crystal  stability  and 
pressure-induced  phase  transformations  in  solids.  These  achievements  are 
Important  because  they  establish  a  new  approach  for  predicting,  from 
first-principles,  the  structural  properties  of  more  complicated  systems  where 
theoretical  understanding  has  been  limited  to  empirical  methods  and  experimental 
data  may  be  lacking.  However,  the  limitations  of  the  LSDA  for  magnetic 


materials,  such  as  iron,  are  still  open  to  question7*7In  this  letter  we  assess 
the  accuracy  of  the  LSDA  in  describing  the  magnetic  and  structural  ordering  in 
3d  transition  metals,  based  on  total  energy  calculations  for  fee  and  bcc  iron  as 


a  function  of  volume  and  type  of  magnetic  order. 


A  comprehensive  study  of  the  ground  state  properties  of  metals.  Including  3d 
ferromagnetic  metals,  has  been  reported  by  Moruzzl,  Janak,  and  Williams1  using 
Korringa-Kohn-Ros coker  (KKR)  method  with  a  LSDA  exchange-correlation  potential. 
Rubier2  subsequently  used  the  aug me nted-spherica  1-wave  (ASW)  method  to  calculate 
the  total  energies  of  non-magnetic  (NM),  ferromagnetic  (FM),  and  anti¬ 
ferromagnetic  (AFM)  states  for  both  bcc  and  fee  phases  of  iron  for  a  range  of 
atomic  volumes  near  equilibrium.  The  FM  state  was  found  to  have  the  lowest 

energy  in  the  bcc  case  while  the  AFM  state  is  the  lowest  in  the  fee  case. 
Surprisingly,  the  latter  was  calculated  by  Kubler  to  lie  1184°K  lower  in  energy 

than  the  former,  which  is  the  actual  observed  ground  state  of  iron  at  zero 
temperature.  The  muffin-tin  (MT)  approximation,  used  in  the  ASW  calculations, 
was  suspected  to  be  the  most  likely  cause  of  the  discrepancy. 

To  further  investigate  this  point,  we  have  used  a  self-consistent 
linearized  augmented  plane  wave  (LAPW)  method3  which  is  free  from  any  shape 
approximations  to  the  potential  and  charge  density.  The  method  has  been  tested 
and  found  to  work  well  on  a  variety  of  different  systems.  Specifically,  we  have 


studied  the  NM  and  FM  phases  of  both  bcc  and  fee  iron,  as  well  as  an  AFM  fee 
phase  which  consists  of  alternating  layers  of  up  and  down  spins  normal  to  the 
[001]  direction.  What  Is  done  for  the  first  time  here  is  to  obtain  a  more 
accurate  estimate  of  the  crystal  structure  difference  in  the  total  energy  of 
iron  which,  in  the  past,  has  been  limited  by  the  use  of  the  MT  approximation. 

We  have  found  the  fee  phase  more  favorable  among  the  two  NM  states,  which  is 
consistent  with  the  fact  that  Ru,  the  4d  counterpart  of  Fe,  is  hep.  In  the  bcc 
case  our  results  give  the  normal  FM  state  as  3334eK  lower  in  energy,  6.7Z  larger 
in  equilibrium  volume  and  15Z  smaller  in  bulk  modulus  than  the  NM  bcc  state.  In 
the  fee  case,  we  find  two  metastable  FM  states,  with  the  ground  state  being 
either  AFM  or  NM  (the  two  are  degenerate  within  our  numerical  accuracy).  Most 
notably,  the  AFM  or  NM  fee  phase  has  the  lowest  energy  (867°K  lower  in  energy 
than  the  FM  bcc  phase).  We  believe  the  failure  to  predict  the  proper  ground 
state  of  iron  Indicates  a  fundamental  difficulty  in  the  LSDA  for  describing  the 
magnetic  Interaction  in  3d  transition  metals. 

Before  discussing  our  results  further,  we  summarize  the  salient  points 
regarding  the  theoretical  and  numerical  approaches  we  have  used.  The  LAPW 
calculations  closely  follow  the  methods  discussed  in  Ref.  3  with  the  iron  core 
states  treated  self-conslstently  using  the  full  Dirac  equation  for  the  spherical 
part  of  the  potential  and  the  valence  electrons  treated  in  the  send -relativistic 
approximation  which  neglects  spin-orbit  coupling.  The  self-consistent  charge 
densities  were  evaluated  using  a  histogram  sampling  of  40,  60  and  30  special  ( 
points  in  the  IBZ  for  the  bcc,  fee,  and  AFM  fee  phases,  respectively,  with  an 
artificial  temperature  broadening  of  2  mRy.  All  non-spherical  and  non-MT 
corrections  to  the  charge  density  and  potential  were  explicitly  included  in  the 
self -consistent  iterations.  Careful  tests  of  the  sensitivity  of  the  total  energy 
to  basis  set  size  and  to  the  choice  of  the  LAPW  energy  parameters  were 


perforated.  Overall,  the  structure  differences  of  our  total  energies  have  an 
absolute  error  of  less  than  1  mRy . 

For  the  LSDA  exchange-correlation  potential  we  used  the  parametrization  of 
Vosko  et  al.4  This  potential  Interpolates  between  the  well  known  high-  and 
low-density  limits,  using  the  exact  Monte  Carlo  calculations  of  Ceperley  and 
Alder5  at  intermediate  densities  for  both  the  NM  and  completely  polarized 
electron  gases.  For  intermediate  spin  polarization,  however,  the  spin  dependent 
part  of  the  potential  is  scaled  according  to  the  random  phase  approximation. 

The  total  energy  per  atom  calculated  for  the  bcc  and  fee  phases  is  shown  in 
Fig.  1  together  with  the  atomic  moments.  The  solid  lines  are  obtained  from 
least  square  fits  to  the  Murnaghan  equation  of  state5  from  which  one  can 
determine  the  equilibrium  lattice  constants  and  bulk  moduli.  They  are  given  in 
Table  1.  For  FM  bcc  iron,  our  calculated  equilibrium  lattice  constant, 

5.212  a.u. ,  and  bulk  modulus,  2.66  Mbar  are  in  good  agreement  with  recent 
independent  calculations  of  Hathaway  et  al.7  The  corresponding  experimental 
values  are  5.406  a.u.  and  1.68  Mbar,  respectively.  The  discrepancy  between 
theory  and  experiment  is  larger  than  what  is  typically  found  for  nonmagnetic 
materials. 

The  band  theory  of  ferromagnetism  is  now  well  understood.  The  magnetic 
Interaction  leads  to  a  spin  splitting  of  the  d  bands  and  a  filling  of  the  less 
bonding  majority  spin  orbitals  at  the  expense  of  the  more  bonding  minority  spin 
orbitals.  This  causes  the  volume  to  Increase,  which  in  turn  reduces  the 
conpresslon  the  sp  electrons  feel  (i.e.  reduces  their  kinetic  ener^),  and  hence 
the  bulk  modulus.  Ferromagnetism  is  established  if  the  Stoner  criterion  is 
realized.  The  NM  density  of  states  at  the  Fermi  energy  must  be  reasonably  high 
so  that  the  gain  in  the  kinetic  energy  is  more  than  conpensated  for  by  the  loss 
in  the  exchange  energy.  The  above  conditions  are  satisfied  in  the  bcc  phase.  We 
found  ferromagnetism  lowered  the  total  energy  by  3334*K,  accompanied  by  an 
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expansion  in  volume  (6.7Z)  and  a  reduction  in  the  bulk  modulus  (15Z).  This 
energy  difference  should  not  be  confused  vith  the  Curie  temperature  (Tc  -  1041°K) 
which  is  the  temperature  above  which  the  ave rage  magnetization  vanishes.  Over 
this  temperature  range  the  thermodynamics  is  dominated  by  spin-wave-like 
excitations  (changes  in  the  direction  of  the  exchange  fields)  rather  than 
Stoner-like  excitations  (changes  in  the  magnitudes  of  the  exchange  fields).8 
Thus,  substantial  local  moments  persists  even  above  Tc  in  contrast  to  the  NM 
states  that  we  have  considered.  This  is  supported  by  many  experiments  Including 
the  exchange  splitting  close  to  Tc  observed  in  photoemission. 9 

At  normal  pressure,  experiment  indicates  that  the  bcc  phase  of  iron  is 
stable  below  1183°K  (a  phase)  and  above  1663°K,  up  to  the  melting  point  of  1807°K 
(6  phase).  Between  1183°K  and  1663°K,  the  fee  y  phase  is  stable  but  the 
temperature  is  too  high  for  long-range  magnetic  ordering  to  occur.  Therefore, 
the  magnetic  properties  of  the  y~Fe  at  low  temperatures  are  not  well  understood. 

Early  measurements  on  the  fcc-Fe  alloys, 10  and  small  particles  of  fcc-Fe 
precipitates  from  supersaturated  Cu-Fe  solid  solutions,11  led  to  the  conclusion 
that  the  ground  state  of  fcc-iron  is  AFM.  However,  more  recent  measurements  of 
small-particle  precipitates  from  Cu-Au  alloys  indicate  FM.12  In  addition,  films 
of  fee  iron  grown  on  copper  surfaces  may  be  either  FM  or  AFM,  depending  on  the 
crystallographic  surface  on  which  the  growth  occurs.13  Recent  neutron  measure¬ 
ments  on  bulk  y-fe  at  high  temperature  Indicates  the  presence  of  FM  correlations 
and  a  substantial  local  moment.11* 

At  normal  pressure  our  results  shown  in  Fig.  1  indicate  that  the  NM  and  AFM 
phases  of  fcc-Fe  are  essentially  degenerate.  This  near  degeneracy  is  consistent 
with  the  low  Neel  temperature  observed  in  fcc-Fe  alloys,  but  is  in  contrast  to 
Kubler's  calculation  which  placed  AFM  fcc-Fe  1183°K  lower  in  energy  than  NM 
fcc-Fe.  Our  results  are  supported  by  the  fact  that  our  calculations  show  that 
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the  total  energy  of  AFM  fcc-Fe  converges  to  that  of  the  NM  fcc-Fe  when  the  atomic 
moment  vanishes.  Unlike  the  FM-NM  transition  in  bcc-Fe,  there  is  very  little 
change  in  the  equilibrium  volume  and  bulk  modulus  between  AFM  and  NM  states  of 
the  fee  iron.  This  may  be  due  in  part  to  the  small  atomic  moments  (0.64  pB)  in 
the  AFM  fcc-Fe  conpared  to  that  of  the  FM  bcc-Fe  (2.08  pB).  The  small 
differences  in  the  atomic  moments  from  previous  calculations1 »2 ,15,16,  17 
reflect  (1)  different  lattice  constants,  (2)  different  choices  of 
exchange-correlation  potentials,  (3)  semi -relativistic  effects  which  tend  to 
lower  the  sp  states  relative  to  the  d-bands,  and  (4)  the  use  of  the  MT 
approximation  in  some  cases. 

In  good  agreement  with  earlier  calculations,2*16  our  total  energy  as  a 
function  of  atomic  volume  for  FM  fcc-Fe  exhibits  two  local  minima,  with  the 
magnetic  moment  increasing  abruptly  between  them.  Thus,  our  results  support  the 
idea  that  FM  fcc-Fe  can  exist  in  two  states  with  only  slightly  different 
energies:  a  low-volume,  low-spin,  large-bulk  modulus  state,  and  a  high-volume, 
high-spin,  small-bulk  modulus  state.  The  coexistence  of  two  spin  states  in  FM 
fcc-Fe  has  long  been  postulated  to  explain  the  invar  anomalies.  The  origins  of 
the  two  spin  states  has  been  discussed  in  terms  of  features  of  the  density  of 
states,16*16  while  Bagayoka  and  Callaway15  have  identified  the  abrupt  transition 
with  the  occupancy  of  a  nearly  flat  portion  of  the  d-bands  and  the  disappearance 
of  the  corresponding  Fermi  surface. 

Our  results  can  also  resolve  a  small  discrepancy  regarding  the  transition 
from  the  low-spin  state  to  the  high-spin  state:  Kubler2  found  the  transition 
between  r8  «  2.64  and  2.65  a.u.  and  a  moment  of  1  pg  while  Bagayoka  and 
Callaway15  reported  a  larger  value  of  rg  (2.71)  and  moment  (1.52  pg).  We  found 
metastable  solutions  for  both  high-  and  low-spin  states  between  rg  -  2.66  and 
2.68  a.u.  In  the  low  spin  phase,  the  corresponding  atomic  moments  increase 
sharply  from  1.35  to  1 • 56  pg.  The  possibility  of  more  than  one  solution  in  a 


self -consistent  calculation  was  first  discovered  by  Janak19  in  ferromagnetic 
cobalt. 

As  can  be  seen  from  Fig.  1,  the  atomic  moment  decreases  with  decreasing 
volume,  which  Indicates  that  there  is  a  transition  back,  to  the  NM  state  as  the 
volume  shrinks,  the  band  width  increases,  and  the  electrons  become  delocalized. 

On  the  high  volume  side,  there  is  another  transition  in  fcc-Fe  from  AFM  to  the 
high-spin  state  of  FM  near  rg  -  2.69  a.u.  This  is  consistent  with  the 
experimental  observations  that  fcc-Fe  precipitates  in  Cu  (rg  ■  2.67)  are  AFM11 
and  fcc-Fe  precipitates  in  Cu-Au  alloys  (rs  *  2.78)  are  FM.12 

The  major  result  of  the  present  study  is  that  both  NM  and  AFM  fcc-Fe  lie 
lower  in  energy  than  does  FM  bcc-Fe  by  867°K.  This  emphasizes  a  shortcoming  of 
LSDA  in  incorrectly  describing  the  effect  of  the  Hund's  rule  correlations  that 
are  responsible  for  the  formation  of  atomic-like  local  moments.  Within  LSDA,  the 
exchange  and  correlation  potentials  are  determined  locally  based  on  the  results 
of  a  homogeneous  electron  gas  of  the  same  spin  density  at  that  point.  The 
determination  of  these  local  potentials  makes  no  distinction  as  to  whether  the 
spins  of  the  nearest  neighbor  atoms  are  parallel  or  anti-parallel  and,  therefore, 
some  of  the  correlations  between  local  moments  on  different  sites  are  neglected. 
It  is  possible  that  the  3d  electrons  in  iron  are  localized  enough  for  this  effect 
to  become  important.  Recently,  one  of  us20  has  studied  the  effects  of  a 
non-local  exchange  potential  and  found  a  sizeable  correction  to  the  exchange 
splitting  in  nickel.  Similar  attempts  in  iron,  however,  failed  because  the 
correction  is  again  large  in  spite  of  the  good  agreement  between  the  uncorrected 
LSDA  energy  bands  and  experiments.  These  results  may  be  an  indication  of 
substantial  cancellation  between  non-local  exchange  and  non-local  correlation 
corrections,  which  must  be  treated  on  an  equal  footing.  Hopefully,  the  present 
letter  will  stimulate  more  work  along  these  directions. 
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Table  I 


Total  Energy  per  atom  (E)  at  equilibrium  lattice  constant  (a),  bulk  modulus  (B), 
and  atomic  magnetic  moment  (y)  of  NM,  FM,  and  AFM  states  of  bcc  and  fee  Fe. 


E 

a(a.u.  ) 

B(Mbar) 

m(Ujj) 

bcc 

MM 

-2541.173 

5.10 

3.14 

- 

FM 

-2541.194 

5.21 

2.66 

2.08 

fee 

NM 

-2541.200 

6.38 

3.44 

- 

AFM 

-2541.200 

6.38 

3.23 

0.64 

FMX 

-2541.200 

6.37 

3.90 

0.00 

fm2* 

-2541.175 

6.81 

- 

2.47 

*There  is  no  equilibrium  solution  for  the  high-spin  Btate  of  fee  FM2  iron  in  the 
sense  that  there  is  no  volume  for  which  the  pressure  vanishes. 
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